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چکیده هدف از این مقاله تحلیل هم‌زمان یک صفحه‌ی دایره‌ای انعطاف‌پذیر و نشیم نآن که یک محیط نيم‌پینهایت بارفتار ایزوتروپ 
جانبی است, می‌باشد به‌طوری که مجموعه‌ی تحت اثر نیروی قائم متقارن محوری موثر بر صفحه قرا رگرفته‌است. در بررسی تحلیلی جزئیات 
اند رکنش بین صفحه‌ی انعطاف‌پذیر و نی‌فضا در اين مقاله ناحیه‌ی تماس صفحه و یم‌فضا بدون قابلیت تحم لکشش در نظ رگرفته شده‌است. 
به‌منظور حل, با اعمال نیروی حلقوی قائم بر دو محیط نیم‌فضا و صفحه به‌صورت جداگانه, تواب م گرین تغییرمکان قائم این دو محیط در 
دستگاه مختصات استوانه‌ای به‌دست می/ید و با در نظ رگرفتن شرایط پیوستگی تغییرمکان‌ها در ناحیه‌ی تماس, معادلات انتگرالی دوگانه‌ی 
حاکم بر تتش تماسی مسأله به‌دست می‌آید. برای حل این معادلات از روش المان‌های محدود استفاده می‌شود به‌طوری‌که با به‌کارگیری 
المان‌های رینگی‌شکل و با تغییر تدریجی اندازه‌های المان, شرایط یکنواختی و تکینگی مسأله قابل مدل‌سازی می‌باشد. به‌منظور تأیید صحت 
انجام محاسبات در این مطالعه, نتایج به‌دست آمده با مقالات پایه‌ای موجود در دو حالت اند رکنش حطی و غیرنعطی صفحه‌ی دایره‌ای 
انعطاف پذیر با نیمفضا مقایسه شده‌است. در ادامه برحی نتایج جدید به‌دس تآمده در این مطالعه به‌منظور نشان دادن اثرات درجات مختاف 
ناایزوتروپی نیم‌فضا به‌صورت جداول و گراف‌هایی ارائه شده است. نتایج این مقاله برای مدل‌سازی واقعی شالوده‌های مستقر بر خاک با رفتار 
ایزوتروب جانبی استفاده می‌شود. 

واژه‌های کلیدی نیم‌فضای ایزوتروپ جانبی صفحه‌ی دایره‌ای انعطاف‌پذیر, متقارن محوری معادلات انتگرالی دوگانه. روش المان‌های 


زین 


۷۲ ۷ )۳۱۵ رماه عاصاه۲۱ مها متا ما۲۱ اه ممتاهونامع ۱۲ مزا امه 
۲۲۵۱/6 15010۲0016 


60. 001-7 ۱ ۳ 2 ۹ 


٩۱۱۵ 6۲۵6۱06۱0۲ ۵ 6 ۶6‏ ۱۱۵۵۵ 0 که ۷ 5۵۵66-]۱۵ ۵06زا ۳۵۹۵۵۹۵ ۸ ]0و۳ 
۵ ۱5 ۳۱۵/۲5۹۵۵۵۵ 6 0۳0 عماج ۳6 ]0 0۲60 60۳۵۵6۲ 1۳6 .عاهام ۲مابدنه ۱۵۱۵ ۵ کیک 066] ای 
عم ۷6۲۵۱6۵۱ 0 « 0۵6۲60 15 090680۴ ۲۳6 .(۱6۸۵10۴8۱655) ۵۱060 ۵0 کعع۳۱6۲0۴۱ 0 ۲۵ 0] 
۰ 115 ]0 0۲961 ۲۱6 و 5۵5 ۸۵0 ۱656 ]0 ۱۱۵۳۵6۵۵ 6 ۵ و0۳۵ 16۵۱60 .۱0۵۵0 0305۳۳۱۱۱۵۵ 
۱69۵ 0۵ ,۱۵۹۵66 و6۵ 6 0۵ ععما عط] اه ۵ معا 2۲۵۵۳85 ۱۵۵0 ود ]0 کی ۵ ۲۷۱۵۲ 
۵۵ ۲0۱9۸۱۵۲ 600۱۴۵1۵ 6 ۱۸۱۵061 0 0۳۵۵8۵0 ۵۲6 ۱869۵۵185 50۳6 ۷۵۲ 008000۵0۴0۱۵0 ک0 690 
۷۵۲۱65 9۲۵0۵۵0 5۱26 ک ۷۳۱۵۲۵ ,۳۱۵۵۱۵۵0 616۳06۳1۲ ۵6 ۲۱9-6۵۵6 ۱6 ]0 ۱۳۵۵۲۵۵۲۵۵۵۴8 ۳6 ۲۷۷۱۵ ۲۵۳۵۱۵۱۰ 
۱6 0۲ ۷۵۱۱۵۲ ۲۳۱6 5۳۱۵۵۵۱۰ 6۵0 ۲مابودنه ۵و۵ ۲مابوع۲ هط وداسنطاویی ]0 6020916 ۵0۲۵ ۱۷۷۵ 
۲ 6 ۱۷۷۱۱۱ 2026۲ کاا ]0 کآیاعع۲ ۱6 و م6۵ ۵0ج کز ۱۱۵۵0۵0 آمع۵ 00 م0۳80 ۳6 ]0 
۷ 50۱6 .۱۵۱-۹۵۵66 00۱0 عاهام ۲مابمه ]0 ۳۵۲۵۵۵۵۲ ۵۳۵۱۳۵۵۲ 0۳00 ۳6۵۲ ۵ ]0 6۵665 ۵۶۵۳۳۱۵۲ 0 

۰ ۱۵ ]0 0۳0۱500۲0۵ 16 ]0 05۴666 ۲۱6 0۲۲۵ ۱0 ۵۳6۹6۱۵6۵ 0۳6 005 کی ]از 


احتعماص1 اهناند رمتتامصصطدرروتن لقلند رماقاظ تقانمتن ماطا۳۱۵ رمموم۵-]۳۱2۸1 15010010 «اموز6 1۳2۵5۷ 16۱۷۵۲۵ 
۰ ۱۵۱۵۱ م1 رقصمتاهتا۴0 


ک تاریخ دریافت مقاله ٩۳/۲/۱۰‏ و تاریخ پذیرش آن ۹۶/۱/۲۵ می‌باشد. 

(۱) دانشجوی دکتری, مهندسی عمران- سازه دانشگاه صنعتی نوشیروانی بابل. 

(۲) نویسنده‌ی مسئول, استاد. گروه مهندسی عمران, دانشکده‌ی عمران, دانشگاه تهران. 

(۳) استادیار گروه مهندسی عمران. دانشکده‌ی عمران. دانشگاه علوم و فنون مازندران بابل. 


۱۳۶ 


مقد مه 

بسیاری از مصالح در طبیعت و نیز ساخته‌های مصنوعی 
رفتار ایزوتروپ جانبی دارند. از آن‌جمله می‌توان به 
خاکی زير ساختمان‌ها و صفحات چندلایه اشاره کرد. 
اهمیت بررسی پاسخ این مصالح از دیرباز مورد توجه 
بوده است به‌طوری که میشل در سال ۱۹۰۰ میلادی به 
بررسی یک نیم فضای ایزوتروپ جانبی تحت نیروهای 
سطحی دلخواه پرداخته ات [24]. لخنیتسکی در سال 
۰ محیط ایزوتروپ جانبی را در حالت متقارن 
محوری و بدون پیچش در نظر گرفته و معادلات درگیر 
حاکم بر مسئله را با معرفی یک تابع پتانسیل به‌صورت 
مجزا و قابل حل در آورده‌است [22]. هو محیط 
ایزوتروپ جانبی را در حالت کلی مورد توجه قرار 
داده و تابع پتانسیل لخنیتسکی را برای حالت کلی 
تکمیل کرده‌است [21]. این تابع هم‌اکنون در ادبیات 
مکانیک محیط‌های پیوسته با رفتار ایزوتروپ جانبی 
تحت نام تابع لخنیتسکی- هو- نواکی نامیده می‌شود. 
محیط با رفتار ایزوتروپ جانبی توسط دیگران هم‌چون 
الیوت[11], ایوبنکس و استرنبرگ[18] پن و چو [27] 
و ونگ و ونگ [31] نیز در حالت استاتیکی بررسی 
شده‌است. دینگ و همکاران در سال ۲۰۰۳ مطالعه‌ی 
کاملی بر روی رفتار محیط ایزوتروپ جانبی داشته‌اند. 
آنها اثر بارگذاری‌های مختلف از جمله تنش‌های ناشی 
از نیرو و حرارت سه‌بعدی را بر روی این محیط 
بررسی نموده‌اند [10]. این محیط در حالت دینامیکی 
توسط اسکندری قادی [12] و اسکندری قادی و 
همکاران [13] و دیگران مورد توجه قرار گرفته است. 

بسیاری از محققان در زمینه‌ی علوم مهندسی و 
ریاضیات بوده‌است. دامنه‌ی کاربرد این مسائل بسبار 
زیاد است و از آن جمله می‌توان به طراحی فونداسیون 
در سازه‌های شتگیرء مهندسی روسازی. انش هباج 
نعیین میزان تغییرمکان در سطح تماس دو محیط 
اندرکنش بین دو سطح و کاربرد آنها در علم مکانیک 
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اشاره کرد. بسته به میزان سختی نسبی مصالح در سطح 
تماس آنهاء هندسه‌ی ناحیه‌ی تماس و بارگذاری وارد. 
مسائل غیرخطی بسیاری می‌تواند مطرح شود و شامل 
تماس سطوح صلب- الاستیک و الاستیک- الاستیک 
می‌باشد که توزیع تنش تماسی در آنها ممکن است 
به‌صورت یکنواخت يا به‌صورت سینگولار باشد. در 
زمینه اندرکنش صفحه‌ی صلب و محیط ایزوتروپ 
جانبی در حالت استاتیکی و دینامیکی تحقیقاتی 
صورت گرفته‌است که می‌توان از تحقیقات اسکندری 
قادی و اردشیر بهرستاقی [۲] و [17 ,16 ,15 .14] 
به‌عنوان نمونه نام برد. در یک حالت خاص از این گونه 
مسائل می‌توان به بررسی مسئله‌ی تنش تماسی در 
سطح بدون قابلیت تحمل کشش بین یک محیط 
نک کی تک همطاف ای از ناد 
متمرکز مرکزی قرار گرفته بر روی آن پرداخت. 

تحلیل یک صفحه‌ی دایره‌ای و فونداسیون وینکلر 
(۱۷10110) توسط تیموشنکو و ووینوسکی در سال 
۹ توضیح داده شده‌است [30]. ویتسمن در سال 
۹ یک راه حل تقریبی برای تعیین شعاع تماس بین 
یک صفحه‌ی الاستیک و یک نیم‌فضای الاستیک اراشه 
دام [ 92 هه بخطا لت ره هدن مایا ور 
سطح تماس بدون قید و تحت اثر یک بار متمرکز در 
وسط در نظر گرفت. هم‌چنین او در سال ۱۹۷۱ به 
بررسی یک تير اویلر- برنولی (۳016-8600[1) تحت 
بار متمرکز متحرک قرار گرفته بر روی یک فونداسیون 
توت فاناش هل کشتی برواست ۳ و راو ای 
مطالعه به اين نتیجه رسید که میزان دامنه‌ی بارگذاری و 
سرعت حرکت بار بر روی ناحیه‌ی جدایی تیر از 
فونداسیون تأثیر بسزایی دارد. او در ادامه‌ی مطالعات 
خود. در سال ۱۹۷۲ به محاسبه‌ی تقریبی فاصله‌ی 
نقطه‌ی اثر بارگذاری روی یک تیر الاستیک قرارگرفته 
بر روی یک نیم‌فضای الاستیک بدون خمش پرداخت 
[34]. در ادامه گلدول و اییر در سال ۱۹۷۶ به بررسی 
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دایره‌ای الاستیک و نیم‌فضای ایزوتروپ پرداختند [19]. 
آنها بارگذاری وارد بر روی صفحه را به‌صورت 
گسترده و متقارن در نظر گرفتند و برای تحلیل نیم‌فضا 
از تبدیل انتگرالی هنکل و برای تحلیل صفحه از تثوری 
کریشهف (0:0108 و ریسنر (615900۴) استفاده 
کردند. در ادامه آنها با استفاده از پیوستگی تغییرمکان 
در نقاط مختلفی از ناحیه‌ی تماس. مقدار تنش تماسی 
را برحسب بارگذاری وارد و سختی نسبی بین صفحه 
و فونداسیون تعیین کردند. سیلیپ [6 ,4 ,3 سیلیپ و 
همکاران [5] سیلیپ و دمیر [8 ,7] و سیلیپ و گولر 
[ف] طایح رو انا استا تن رورا ز عاشتات 
هارمونیک دینامیکی بر روی تي صفحات دایره‌ای و 
مستطیلی قرار گرفته بر روی فونداسیون‌های بدون 
قابلیت تحمل کشش از نوع وینکلر و الاستیک انجام 
داده‌اند. به‌عنوان مثال» سیلیپ و گولر در سال ۲۰۰۷ به 
تعیین پاسخ استاتیکی و دینامیکی صفحه‌ی دایره‌ای 
قرار گرفته بر روی فونداسیون بدون اصطکاک با 
تعریف دو پارامتر پرداختند. آنها در تحقیقاتشان» تنش 
تماسی ۲ را به‌صورت دو پارامتری به شکل 
۷-۷ تعریف کردند که ,5 و ع دو 
پارامتر برای مدل‌سازی اثر خاک زیرین فونداسیون و 
تغییرمکان فونداسیون می‌باشد. پک و همکاران در 
ال ۳۰۸۸ به خی کامان اش تون فالخ قفا 
کشش یک صفحه‌ی انعطاف‌پذیر دایره‌ای قرار گرفته 
بر روی نیم‌فضای ایزوتروپ پرداختند. آنها مسأله‌ی 
مایت تا فد این ی و تور ماو لاش 
انتگرالی و شرایط مرزی و پیوستگی حاکم بر مسأله‌ی 
مدل‌سازی کردند و شعاع ناحیه‌ی تماس را به‌عنوان 
یک مجهول اضافه به روش سعی و خطامحاسبه 
نمودند [26] 

تا جایی که نویسندگان مطلع هستند. اندرکنش 
صفحه‌ی انعطاف‌پذیر و محیط ایزوتروپ جانبی بررسی 


تستالته اشتت: لذ| در این مقاله این موضوع مورد توجه 
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قرار گرفته‌است. بدین منظور ابتدا با اعمال نیروی قائم 
تگن بو بش لیف تیفآنیان سا گوریر 
صفحه تابع گرین تغییرمکان قائم صفحه‌ی دایره‌ای 
به‌دست می‌آید. در ادامه‌ی معادلات حاکم بر نیم‌فضای 
ایزوتروپ جانبی در سیستم مختصات استوانه‌ای شامل 
معادلات تعادل. روابط تنش- کرنش یا معادلات 
رفتاری و روابط کرنش- تغییرمکان بیان شده و در 
ادامه معادلات تعادل بر حسب مژلفه‌های بردار تغییر 
مکان به‌دست می‌آیند. این معادلات شامل دو معادله‌ی 
دیفرانسیل درگیر با مشتقات جزئی می‌باشند که برای 
مجزاسازی آنها از تابع پتانسیل لخنیتسکی استفاده 
می‌شود. در ادامه به‌کمک تبدیل اند رالین هنکل. 
مولفه‌های بردار تغییرمکان در فضای تبدیل‌یافته برای 
نیروی قائم با توزیع رینگی به‌دست می‌آید و با اعمال 
تبدیل معکوس هنکل. توابع گرین تغییرمکان قائم و 
افقی نیم‌فضا به‌دست می‌آید و با کمک تابع گرین 
تغییرمکان قائم صفحه‌ی انعط‌اف‌پذیر و نیز شرایط 
پیوستگی موجود در سطح تماس آنها طریقه‌ی به‌دست 
آوردن معادلات انتگرالی دوگانه مورد بحث قرار 
می‌گیرد. با توجه به انعطاف‌پذیر بودن صفحه ناحیه‌ای 
از صفحه از نیم‌فضای زیرین وابسته به نوع بارگذاری و 
سختی نسبی صفحه و نیم‌فضا جدا می‌شود. شعاع این 
ناحیه به روش سعی و خطا تعیین می‌شود. سپس 
معادله‌ی انتگرالی به‌دست آمده با استفاده از روش 
المان‌های محدود و به‌کارگیری المان‌های رینگی با 
توزیع ثابت فشار حل می‌شود. از نتایم حاصل فشار 
زیرصفحه‌ی انعطاف پذیر تعیین شده و از آن اندازه 
نیروی تماسی در اثر تحریک وارد به صفحه به‌دست 
می‌آید. نشان داده می‌شود که نتایج به‌دست آمده از این 
روش برای مواد ایزوتروپ در حالتی که سختی صفحه 
به سمت بی‌نهایت میل کند (صفحه‌ی صلب) بر نتایج 
قبلی ارائه شده توسط لوکو و میتا [23] و در حالتی که 
سختی صفحه محدود است (صفحه‌ی انعطاف‌پلذیر). 
بر نتایج ارائه شده توسط پک و همکاران [26] منطبق 
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نشان داده می‌شوند. 


چگونگی حل سأله 
صفحه‌ی حلقه‌ای انعطاف‌پذیر به شعاع داخلی ,6 و 
شعاع خارجی 2 مستقر بر نیم‌فضای ایزوتروپ جانبی 
تحت اثر بارگذاری قائم متقارن محوری ٩,(‏ و تنش 
تماسی (0,)۲ در نظر گرفته می‌شود (شکل ۱ و۲). ,ث 
مدول یانگ ,۷ ضریب پواسون و ضخامت صفحه 
می‌باشد. ,0 شعاعی است که تا قبل از آن صفحه و 
نیم‌فضا با هم در تماس بودند و بعد از آن صفحه از 
نیم‌فضا جدا خواهد شد. فرض می‌شود که تثوری 
کریشهف بر صفحه حاکم و هم‌چنین سطح تماس 
صفحه و نیم‌فضا بدون قابلیت تحمل کشش باشد. در 
ابتدا معادله‌ی حاکم به‌صورت کلی برای صفحه‌ی 
حلقه‌ای تحت اثر بار رینگی اعمال شده در شعاع "۲ 
حل می‌شود و پس از تعیین ضرایب ابت حاصل از 
انتگرال گیری» تابع گرین حاکم بر تغییرمکان قائم 
صفحه‌ی حلقه‌ای به‌دست می‌آید. در ادامه برای داشتن 
این تابع در حالت صفحه‌ی دایره‌ای کامل با میل دادن 
شعاع داخلی مه به سمت صفر ضرایب برای صفحه‌ی 
کامل به‌دست می‌آید و از آن تابع گرین تغییرمکان قائم 
صفحه‌ی دایره‌ای به‌دست می‌آید. سپس معادلات تعادل 
حاکم بر نیم‌فضای ایزوتروپ جانبی در نظر گرفته 
می‌شود و با اعمال بار رینگی» تابع گرین ملفه‌های 
تغییرمکان و تنش به‌دست می‌آید و در ادامه با استفاده 
از شرایط پیوستگی تغییرمکان قائم. تابع تنش تماسی 
صفحه ((),0) با حل معادلات انتگرالی حاکم بر 
مسأله و با استفاده از روش المان‌های رینگی تعیین 
خواهند شد. در ادامه با داشتن تابع تنش تماسی مقدار 
شعاع تماس ,0 به عنوان مجهول اضافه با سعی و خطا 
قابل محاسبه خواهد بود و از روابط تنش-تابع پتانسیل 


و تغییرمکان- تابع پتانسیل, مقادیر تغییرمکان و تنش 


نشربه‌ی مهندسی عمرآن فردوسی 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف پذیر .. 


در نقاط مختلف صفحه و نیم‌فضا به‌دست ی آیقا 


٩ 6,۵,2‏ نیم فضای ایزوتروپ جانبی| 


15, 1 0, 5, 7 
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شکل ۱ صفحه‌ی حلقه‌ای انعطاف‌پذیر مستقر بر نیم‌فضا 


(2 ,۵ ,6 ۵ نیم فضای ایزوتروپ جانبی| 
0 ,0 ۶ 5 1 


۱2 


شکل ۲ اثر تنش تماسی صفحه‌ی حلقه‌ای بر روی نیم‌فضا 


تعیین تابع گرین بار رینگی حاکم بر صفحه 
صفحه‌ی حلقه‌ای تحت اثر بارگذاری قائم متقارن 
محوری ٩,)2(‏ و واکنش بستر (0,)2 در نظر گرفته 
می‌شود. معادله‌ی دیفرانسیل حاکم بر صفحه برای 
تغییرمکان عمودی (),۱۷با فرض تغییرمکان‌های 
کوچک در دستگاه مختصات استوانه‌ای به‌صورت زیر 
می‌باشد [29]: 

83 1 0/۷ 1 ۲ ( 


ی 


۳۹ 
تِ 
)0 


سال بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


قاس مگرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


که در آن ۲ 
۳8 ۳ 
90 پتس- 
( ۷ - 1201 
و (),0-(۵0(<0,6 برآیند بارگذاری قائم وارد بر 
هر نقطه از صفحه می‌باشد. با در نظر گرفتن قرارداد 
برشی شعاعی صفحه به‌صورت زير می‌باشند [29]: 


1۳, ۰ 0۳, )( 


۱۳۷ 


تمتظورن نهدسست آوردق تابع گرین صفحه به‌حای 
بار رینگ دی با شدت 
0۵0 بار رینگی عمودی ب ۳ 


به‌صورت زیر و مطابق شکل (4) در نظر گرفته می‌شود: 
6 ۵0-1 


270۳ 


در شعاع ۳ 


< (0)۳ 
که در آن 8 تابع دلتای دیراک می‌باشد. 

با تقسیم دامنه‌ی شعاعی صفحه به دو ناحیه با 
فصل مشترک در شعاع به‌صورت زیر 


(۲۳ >> وه 22 > 60 0 |(0,)) < ,3 


۳ (- - ۷ 
۳( 7 ۳ ۲ ۹9 )0( 
۹2 ۵۱ >> 2,۲ >0 > 0 |(2,0)) > و 
7 یت ٩‏ نس اف ۳۳ - ۵) ۷ 
۴ 1۳۰ ۳ 69 0 
90 ۷ 
۲ ۱۲۵۱۲ چ س 
۴ 0 
۱ ۲ | مسلط 
(0:)۲ 
واه 2 ۲ 


شکل ۳ المان دیفرانسیلی صفحه‌ی تحت اثر بارگذاری قائم وارد بر آن در مختصات شعاعی 


22 7 


/ 7 ۳ ی 


نیم فضای ایزوتروپ جانبی 


۳ ِ 


/ / ش 3 طً 


شکل ۶ صفحه‌ی حلقه‌ای انعطاف‌پذیر تحت اثر بار رینگی وارد بر آن در شعاع "1 


سال بیست و هفتم» شماره‌ی یک» ۱۳۹۶ 


نشریه‌ی مهندسی عمران فردوسی 


۱۳۸ 


جواب معادله‌ی دیفرانسیل مرتبه‌ی چهارم حاکم 
بر صفحه در این دو ناحبه به‌صورت زیر به‌دست 
می‌آید: 
( . رده می6+ تب +تصاین+تطا < ( ,۷ 
( بت می6+ تب +تصین+تطا کی (0 ,۷ 
در روابط فوق مقادیر ,) الی ,نک ضرایب ثابت 
تخاضل. از انیگرال کیزیف باشکن که با اتتفاده اززقن ایط 
مرزی و پیوستگی حاکم بر صفحه تعیین خواهند شد. 
با توجه به توضیحات فوق و با در نظر گرفتن 
شرایط گشتاور خمشی و نیروی برشی صفر در ۵ < 2 
تغییرمکان و گشستاور خمشی صفر در ,۲-8 و 
پیوستگی تغیبرمکان. شیب. گشتاور خمشی و نیروی 
برشی در ۲۲۳ شرایط مرزی و پیوستگی حاکم بر 
صفحه‌ی حلقه‌ای به‌صورت : 
,0 2 (۵) ۷1۶ ,0 < (مه) !۷1 ,0 < (۵), ۷ 
( )۷۶ 2 ( ۷۲ ,(۱)۲ ۷ 2 ( ۷۱6۲ ,0 < (۵) ۷۶ 
8 


رت 


(0*ب 4 - 


۷۱۱-۷6 
1۹ 


7/9 
)۱۰ 
اکتا با جای گذاری روابط (۰۳ (ع (۸) و )۹( در 
روابط فوق دستگاه معادلات حاکم بر ضرایب ب) الی 


1۳-27 اف 


) به صورت زیر نوشته می‌شوند: 


,20 )+ توبن + مدط ین + مفطا ره 6 


1- ,۷ 
ِِِ_ و +((,۷ +1),ه صا2 +3+ ,)6 
0 
,2 (,20,)1+۷ 
1 ۷ 
اب ۵ +((,۷ + صا2 + 3+ ,یه 


طر46 


2 (,۷ +0 ,20 
0 18۲6-01۲0 ۲۳ 
0-0 +( 0-:0) ۲2 
+(3+ +( 0+۷ 210۲( -6) 
1 ۷ 
و و( و) 
اس 

٩ +2۳6 -0,( 0,‏ - د راما 0۳ -6) 

1 


4 
م0 < - (ر)- و)- 
271 1 


(۱۱) 


نشریه‌ی مهندسی عمرآن فردوسی 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پدیر ... 


با حل دستگاه معادلات فوق. ضرایب ً) الی 


ول) برای صفحه‌ی حلقه‌ای به‌صورت: 


1 


و 
عمط تطا لا 6 ]مق 
۷ 0( مه 1910 
5 و2 1 
۵ ره ۱0۳ 206 | (۱+,(- 16۳12۵7 ِ 
| () یب 2۵6+ (۷۱+3) و3 
2 2 ك 
+(+ ,2+ ,1۳ م2 - 4 
مقط 2 وق ولا )م23 
+((ععاجل ۷ -1)+(1- مققا 
+ (لبم- ,۷+( 31-۷ 2] ‌ 
|(3- 2۷+ 0۳ 
6-0 
۱ _ (0- + (10- تص)(1+ ۳( 2۸0) 5 
2 0 رت ۹۳0 ۱ 
۳ 7 1 
+ رفص ص1۳۷ و2 | 71۳ ۲۳ 
|( ه) 
2 2/2 1 
+ رده 9 5 (1- ,6( مه 1671۵ ی 


2 2 / 

2, 1۳2-۷, + 

+ (مه20- + ,2۷+ 0۳ عطا مفصا2 

4 + تلا 210 1۳۲ 

‌ (( شرب ( ,0-7 ۳ ۳ 

+ (3- 2۷ ۷ ها 

-ِ (تیس 0۲1 ,21) ۳ 

[(9+ ره 

(۱۲) 

به‌دست ات ات با جای‌ گذاری این ضرایب در 
حلقه‌ای در شعاع ۲ ناشی از بار رینگی فائم در شعاع 


۳ به‌صورت زیر به‌دست ری نگ 


ستیالج بیست و هفتم» شماره‌ی یک. ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


(0- 116( بت + ۲7 ب6 +1۲ رن +191 ۲۳ ۳ 
116-۰( + ۳ ب +طتصآ م0 + ها ۲۳ )+ 
۱۳( 

که در آنْ ()11 تابع هویساید است و ضرایب 
ب) الی بن) در رابطه‌ی (۱۲) تعریف شده‌اند. اگر 
شعاع داخلی و2 به‌صورت حدی برابر صفر در نظر 
گرفته شود. ضرایب ,) الی ون) برای صفحه‌ی 
دایره‌ای به‌صورت زیر به‌دست مش ام 


[- 
,20 رن ی 

9 

(,1-۷) ۳۳ +(1+ ربج (ا جر تها) 227 
,20 هناد 7 < پل 

(+ ,0۷ 167122 
۰-1 ۷) ۳7 1 
۱ / 2 رل ,سدع 6 6-0 
(+ ,0۷ 167122 911 

ِِ< 8 ( 


581 
۹13 
با جای‌گذاری ضرایب فوق در روابط (۸) و )٩(‏ 
تابع گرین تغییرمکان قائم صفحه‌ی دایره‌ای در شعاع ۲ 
ناشی از بار رینگی قائم در شعاع "۲ به‌صورت زیر 
به‌دست ره 


2 ۱ 2 
و(« ع+ + +( عصا) دوم 


2 ۷" 
,)167۳۵ 
2 
ی ۲۱ 
زین 
۱ 2 1 2 2 
11)۲-۲(۰ >« ۹ ۳ / ,0 سم _1 
8 م۱6 ناه 
(۱۵ 


با استفاده از تابع گرین تغییرمکان قائم صفحه 
دایره‌ای (رابطه‌ی ۱۵) و نظریه‌ی تقابل بتی [20] 
تغییرمکان عمودی این صفحه تحت اثر برآیند نیروهای 
ناشی از بار قائم خارجی و تابع تنش تماسی به‌صورت 
رابطه‌ی انتگرالی زیر بیان می‌شود: 


۷, )( < ۷, )0( + ۱۳ 6:1 (0, )( - 0, (۰ 


0-1 <>2 
(۱1) 


سال بیست و هفتم» شماره‌ی یک ۱۳۹۶ 


۱۳۹ 


تعیین تابع گرین بار رینگی حاکم بر نیم‌فضا 

نیم فضای ایزوتروپ جانبی که محل نشیمن صفحه 
می‌باشد طوری در نظر گرفته می‌شود که محور 
ایزوتروپی آن عمود بر سطح نشیمن صفحه باشد. 
دستگاه مختصات استوانه‌ای (0:۲,0,2) طوری قرار 
داده می‌شود که در آن امتداد 2 موازی محور ایزوترویی 
باشد و نیم‌فضا با 0< 7 تعریف شود. معادلات تعادل 
برحسب تغییرمکان در شرایط متقارن محوری و برای 
حالتی که نیروهای حجمی وجود نداشته باشد 
به‌صورت زیر نوشته می‌شود [122: 


67۲ 100 00 


0 
0 )+ توح ( 2 1 تچار۸ 


2 2 
۱ ار ۳۹ 


:2۳۵۹ مرح ۲ 
20 22 

کت ات زو ۱۰۱۳ 
تج مج( دش 


(۷ 


در اين رابطه ۷ و ۷ مولفه‌های بردار تغییرمکان 
به‌ترتیب در امتدادهای ۲ و 7 هستند و ۸ نت ایت 
ارتجاعی می‌باشند که تانسور تنش وه را به تانسور 
کرنش وه مربوط می‌سازند [22: 


قورظ ‏ معیرظ ع و6 میرقور گم رطع 6 


223 


قییظ + قورع ره 


و2 م0 ور 
)۱۸( 
کی فیس بو رت 2 رتاش کر 
تعریف می‌شود. این ضریب در ادامه به‌منظور نرمال 
کل یایاده شوه فا صراتب: ارشتاعی 


تون |22 
م۳۲ 
و و 11 
1-0-07 (1۳ ۵-22( +۱) 
3 ۳0-۷۰ 
لا < 20 ۳ بل < پر ‌ ‌ < وی 
(۵ +0 1-0-2۳ 
۱۹( 


نشریه‌ی مهندسی عمرآن فردوسی 


در رابطه‌ی (۱۹) ۳ معرف مدول یانگ در 
صفحه‌ی ایزوتروپی؛ ۳ معرف مدول یانگ عمود بر 
صفحه‌ی ایزوتروپی. ۷ معرف ضریب پواسون در 
صفحه‌ی ایزوتروپی. 1 معرف ضریب پواسون امتداد 


عمود بر صفحه‌ی ایزوتروپی نسبت به هر امتداد در 


صفحه‌ی ایزوتروپی می‌باشد. ۲ معرف مدول برشی در 
صفحه‌ی ایزوتروپسی و ۲ معرف مدول برشی در 
مواد همسان (ایزوتروپ) ضرایب ارتجاعی ۸ 
رت ور 
2 -2001 
0 کب 
1-0 1-۷0 
< م۸ ع بیش 


< بیظ < رر۸ 
۳۲۰( 

در روابط فوق. ۸ معرف مدول برشی و نا 
معرف ضریب پواسون می‌باشد. 

الا اوا ۱۸ ام ی وت 
دفراتفین هرگییا فان جر بسن تور 
مجزاسازی این معادلات از تابع پتانسیل اسکالر 
لختسکی (ظ) ععلت کافال بزدن فالخ ارن 
محوری و سادگی استفاده می‌شود. مژلفه‌های بردار 
تغفییرمکان برحسب تسایع پتانسیل ۲ در دستگاه 
متعضات انطرانه‌اق و دو خالت انسهایکین بصضووت 


زیر نوشته می‌شوند [31]: 


2 2 93 
(۲۱) .۱۳ 0+ ۷ (1+0)< 7 , 0- < 1 
رح ۱ ( 202" 
که در آن: 
9 هس ۰ 1 5 
معط مه مه 
2 
2 ۱ مم 8 
۲ 0۳۳ 0 +1 
(۲۲( 


با قرار دادن روابط (۲۱) در معادلات تعادل (۱۷) 


زیر در می‌آید: 
22 
690 ۰ < (۳)۲,2 ۷ ۷ 


نشریه‌ی مهندسی عمرآن فردوسی 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف پذدیر 5 


که در آن: 
شن (2با<) , 1 7 
027 3 


و پارامترهای ۹ و و8 ریشه‌های معادله‌ی زیر هستند: 
ی ی وبشر ,۸ - پیش و2۸ + جه) ی "شوش 
20 برشررظ 
(۲۵( 

با توجه به مثبت بودن انرژی کرنشی ,5 و 5 
می‌توانند اعداد مختلط باشند اما نمی‌توانند اعداد 
موهومی خالص باشند [22]. به هر حال قسمت حقیقی 
(5 و رو مثبت هستند. با توحه به هندسه‌ی مسئله و 


شرایط مسئله در ۶های بزرگ از تبدیل هنکل مرتبه‌ی 


صفر نسبت به امتداد شعاعی به‌شرح زیر استفاده 


۳0( (ع) !۳2۲ و )۳۳ 
به‌طوری‌که تبدیل معکوس هنکل آن عبارتند از: 
2 (ع)ا ول (ع) ۳ 1 (2 ,۳ 


و ول تابع بسل نوع اول از مرتبه‌ی صفر می‌باشد. 
با اعمال عملگر تبدیل هنکل مرتبه صفر به معادله‌ی 
(۲۳) و استفاده از تعریف (۲۱). معادل هی (۲۳) 


به‌صورت زیر درمی‌اید: 


۲72 72 0 

)۳۸( 0 < (2 نت ۳ و۷ ۷ 
که در آن: 

(۳۹ (2 ,1 رت ( ۷۳-۵6۳4 


| 0 

مرتبه‌ی چهارم با ضرایب ثابت می‌باشند. جواب 

معادله‌ی (۲۸) با توجه به شرایط در دوردست 

به‌صورت زیر به‌دست می‌آید: 

۳ ۳ (ع80 + ۳ ع(ج۸ - («یج) ۳٩٩‏ 
مطابق شکل (۵) فرض می‌شود نیروی استاتیکی 

قائم به‌شدت ()؟ وق تاکز شعاع ۳ در 

0< 7 اعمال می‌گردد. این نیرو به‌صورت زیر در امتداد 

6 تعریف می‌شود: 

۳۱( یر 


79 
271 


شنال بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


قاس مگرجی بندپی- مرتضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


/ 2 ,۵ )۰ نیم فضای ایزوتروپ جانبی ۱ 


آ/ ,0 , 1 ۲ 


شکل ۵ نیم‌فضای ایزوتروپ جانبی تحت نیرو با توزیم حلقوی 
موثر بر حلقه‌ای به شعاع "1 


که تبدیل هنکل آن به‌صورت زیر می‌باشد: 
(۳۲( كت ‌ "0 ۶رع) )۳ 


به‌طور واضح اندازه‌ی کل نیروی ,801 برابر 
0 27 
۳ ] 
121 


6-0 ۲2-0 


شرایط مرزی در 720 با توجه به رابطه‌ی (۳۱) 
و شکل (۵) عبارتند از: 
(:۳( ,0 <- 0:۳ < ۲,2) 6 

۰( ۲,1)] < ( 0:1 < 2 ,)بر 

تور انب ۳ انتد از( ابیسسيي: 
همان‌طور که پیشتر نیز استفاده مایت کلیه‌ی 
ملفههای تانسور تنش و بردار تغییرمکان در 
دوردست (0 جب 2) یا (۰0<- ۲) صفر هستند. 

,20 (02 سل ,0 (2,)رزی حصنا 


0ب 22 2 5 00 وب 2 س 
۰ < (۷۷)۲,2 113 
00 وب 22 1 


(۳۹( 
شرایط مرزی داده شده در روابط (۳۶) باید در 
فضای تبدیل‌یافته نوشته شوند تا با جایگزینی معادلهة 
(۳۰) در آنها بتوان توابع (۸ و (08ظ را به‌دست 


آورد. بدین منظور مولفه‌های بردار تغییرمکان و تانسور 


ال بیست و هفتم» شماره‌ی یک ۱۳۹۶ 


۱۶۱ 


مر بو ۲ 0 5 0 
تنتش در فضای هنکل به‌صورت زیر برحسب تابع ۲۳ 
نوشته می‌شوند: 

(0 


۳0( ون < (۶,2) 1 
17 
۳۷( ۳ تس + ت  )1+‏ («2ع) ۷۳ 
7 
+(0 + اویش - وه به]گ ] وج 


۳ 
۳ ک م0 بیط 


(۳۸( 
۳ 0۱ +1) مب بش (به اه < (2 ی 
(۳۹( 

با جای‌گذاری تابع "۳ در روابط (۳۰) الی (۳۹) 
توابع تغییرمکان و تنش در فضای تبدیل‌یافته نوشته 
مولفه‌های تغییرمکان و تنش برای نیروی حلقوی در 
فضای واقعی به‌صورت زیر به‌دست می‌آیند: 
(ج),1| ۳ 30 یه+ *0ه)۵ه ]و | ( جر ,تا 
(6۰) ۱ 

| )9و وج 0۸+ | 1 به+1) - نیت 1۳ 
(ج) ,|( 80+ *ع)ه) 
)۱ 
«((ه+1) اویش - ورشیه ) مج | - (۳ز2ر) بر 5 
یرهبی۸- [ ۶ 6(ع) ۳ رود ۱ ع()9,۸ 
(ج) ,| ( 90۳ یود ۸/0۳ 
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۱ 
۱ 
۱ 


+ مب هرق ه )مه - یم ]۸ - ( جرج 


وه (ج) ,| ( 130 + ۷ ع/ه) ره +) 
(۳) 


با استفاده از شرایط مرزی (۳۶) و روابط (۳۷) و 
(۲۸) ضرایب ()۸و (8 به‌صورت زیر برحسب 


( )2۳ نوشته می‌شوند: 


نشریه‌ی مهندسی عمرآن فردوسی 


(6) ,۱ (1مج) ۳ 
( 2 - (ع)۸ 
بب۸ ()1 ۳ 
(ج) ,1۱( )1۳ 
(: ما 2 - - (ع)۳ 
و۸ ()1 1 
که در آن: 


,( وک + توفمه - ره +1) < (ج) ,1۱ 
وی( رفیه + "که - به +1) < (ج) ر[1 
م٩‏ < (و) 2 « وه < () ,۸ 

«(به +01 - )یه < (8), 

«(به +1) ی - (ج) میه < (ج) یه 


۸ 

رو( خلت ول ور - 0 +1) < (6) ,۷ 
بب۸ 
۸ 

,۹( 3 0 5 0 +1) < (6) ۷ 
یب 


1)( < 112 )5(۷, )( - 11, )5(۷2)5( 
1) 

با استفاده از توابع گرین مژلفه‌های تغییرمکان و 
تنش نیم‌فضا (روابط ۶۰ الی 4۳) تغییرمک‌ان‌هاو 
تنش‌های نیم‌فضا تحت اثر تنش تماسی صفحه در 
فضای واقعی به‌صورت روابط انتگرالی زیر بیان 


می‌شوند: 
(۷:) 1 (0,2) ,رل 
0 

(6۸) رس (1,2) ,۱۷ 
0 

)۹ (ع رم( برچ أ- (7,) ر6 
0 

(۵۰) ۰( ۲ ,0( 27 و2) مره 
0 


معادلات انتگرالی و شرایط پیوستگی حاکم بر مسأله 
برای در نطر گرفتن شرایط فیزیکی مسأله. فرض 
می‌شود که ناحیه‌ی تماس صفحه و نیم‌فضا به‌صورت 
یکنواخت بدون قابلیت تحمل کشش باشد. هم‌چنین 
به‌دلیل داشتن لبه‌ی آزاد در صفحه‌ی دایره‌ای (شعاع 
0 یک شعاع بلندشدگی (0) وجود خواهد داشت. 


نشریه‌ی مهندسی عمرآن فردوسی 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پذیر .. 


از آن‌جایی که ناحیه‌ی تماس قابلیت انتقال کشش ندارد 
و با توجه به شرایط پیوستگی جابه‌جایی عمودی 
صفحه و نیم‌فضاء روابط زیر در ناحیه تماس برای 
صفحه دایره‌ای حاکم می‌باشند: 
(۵۱ 0 > 1۲ > 0 ,(۲)م ۷۷ < (۲,0) م۷ 


,)۲( < 0, 


برای اطمینان از عدم نفوذ صفحه در داخل 


0>۲ > ۵۲( 


نیم‌فضا در ناحیه‌ی بلندشدگی (2 > ») رابطه‌ی زیر 

نیز باید برقرار باشد: 

(۵۳ ۶۵>> ,0 2 (0,)۲ ,(۲)م ۱ < (۲,0) ,۷ 
اگر تماس صفحه و نیم‌فضا به‌صورت کامل برقرار 

باشد به‌عبارت دیگر اگر کلیه‌ی نیروهای تماسی در 

همه جای ناحیه‌ی اتصال فشاری باشد. 2 در روابط 


(۱) الی (۵۳) جایگزین 0 می‌شود. 
۵ 


به‌صورت زیر نوشته می‌شود: 
اس ( ره( 9 |( رورم ۳ | +۳0 
0 0 
0 >> 0 )0( ۳0:۳ | 
0 
( 
که در آن 1,۳۳ تابع گرین تغییرمکان قائم 
صفحه‌ی حلقه‌ای می‌باشد. رابطه‌ی (۵۶) می‌تواند 


به‌صورت زیر نیز نوشته شود: 
(0), 0-۱۷( ره یمک + رم 2 ]| 
0 


> 0۲ ,۲ ( ۵۹,۲( ۲:۲) ,۱۳ ت- 
0 


(501( 
برای در نظر گرفتن شرایط تعادل کلی صفحه در 
امتداد قائم باید مجموع نیروهای خارجی در تعادل با 
تنش تماسی صفحه باشد. بنابراین رابطه‌ی زیر نیز باید 

برقرار باشد: 
69 2 روت 2 - ۳ 
0 0 


شنال بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


همراه با محدودیت های ذکر شده در روابط (۵۲) 
و (۵۳) معادلات انتگرالی (۵۵) و (۵7) با متغیرهای 
۳,0 و (۷,6۵ تشسکیل یسک سیستم معصادلات 
۱9 حاکم بر مسائل تماسی در حالت متقارن 
محوری را می‌دهد. هم‌چنین برای مسائل تماسی بدون 
قابلیت تحمل کشش مقدار شعاع تماس 0 در نقطه‌ای 
تحلیل مسأله تعیین شود. 
است و به همین علت تحقیقات چندانی در این زمینه 
وجود ندارد. در ادامه با استفاده از روش المان‌های 
رینگی. مقدار واکنش بستر (0,)۳ که مجهول اصلی 
این معادلات می‌باشد به‌دست می‌آید و از آن می‌توان با 
استفاده از روابط مورحود. مقادیر تغییرمکان و تنش 
صفحه و نیم‌فضا را در هر نقطه‌ی دلخواه به‌دست آورد. 


روش المان‌های رینگی 
برای تحلیل نیم‌فضای ایزوتروپ جانبی تحت اثر 
صفحه‌ی انعطاف پذیر دایره‌ای به شعاع 8 نیم فضای 
مورد نظر تحت آ تابع فشار رینگی با شدت ابت 
(1,2,3,....2 < ) ,ره مصوثر بر حلقه‌ای بسه شعاع 
مرکزی 1 و عرض نا قرار می‌گیرد به‌طوری که لا 
تخلاف حعههاین پاش که فست فان باتوی راز 


دارند (شکل‌های 1). به‌منظور انجام محاسبات دقیق‌تر 


۲ 


۱۳ 


بهحصوص در لبه‌های خارجی صفحه عرض هر المان 
از رابطه‌ی زیر محاسبه می‌گردد: 


ما + اب و 
۸ 1 


2 
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(۵۷ . ((,.,12,3ع , 


که در آن 10 ثابت دلخواهی می باشد. به‌ازای 
9-0 تمام المان‌ها عرض ثابتی خواهند داشت و با 
افزایش مقدار ع عرض المان‌ها با نزدیک شدن به 
لبه‌ی صفحه کاهش خواهد یافت. این ویژگی شرایط 
تکینگی به‌حصوص در لبه‌ی خارجی صفحه‌ی صلب 
دایره‌ای را ارضا خواهد کرد. 

در این صورت. تغییرمکان هر نقطه از محیط با 
مختصات (1۷1)7,72 برابر: 


(5۸) 


می‌باشد. به‌طوری که (),۱۷ تغییرمکان نقطه‌ی ۷ 
به‌علت فشار حلقوی رم مور بر حلقه‌ای به شعاع 1 و 
به پهنای با است. از آن‌جایی که ,0 ثابت می‌باشد. 
(۷۷,)۷1 با استفاده از رابطه‌ی (4۸) برابر است با: 


.10 
۷,2 | ۵ - ۷۰۵۵ 0 - (۷) :۱ 
۳3 
1 1 
 12,3,...,2((‏ دشر ۱۳ ۳ > ,19 
از احز 
64٩(‏ 
۳۰ 
8 
0 
1 
۴۳:۳۵ 
با<...جا<با 
۸ ۰ 3 2 4 
۰ 0 ۰ 
۳ 


شکل 1 نیم‌فضای ایزوتروپ جانبی تحت فشار با توزیع حلقوی موثر بر لا[ حلقه با پهنای متغیر در صفحه‌ی دایره‌ای 


سا بیست و هفتم» شماره‌ی یک ۱۳۹۶ 


نشریه‌ی مهندسی عمرآن فردوسی 


۱ 


۱ حح< 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پذدیر .. 


۳4 


شکل۷ نمايش شماتیک فنر معادل نیم‌فضا برای صفحه‌ی صلب دایره‌ای به شعاع 2۸ 


از آن‌جایی که صفحه‌ی دایره‌ای مورد مطالعه 
انعطاف‌پذیر می‌باشد. تغییرمکان نقاط مختلف نیم‌فضا 
در محل تماس با صفحه. غیریکنواخت است و برابر 
تغییرمکان صفحه می‌باشد. با مساوی قرار دادن 
تغییرمککان آ نقطه (0-< 21,2 1۷:۲ ببس 
(1,2,3,...,۷ < ذ) ,() ,۰۷ لا معادله تشکیل خواهد 
شد. به‌طوری که باحل آنها ( مجهول 
(,...,1,2,3 < ) ,.ط تعیین می‌شوند. با در اختیار 
داشتن ,0هاء توزیع تغییرمکان و تنتش در هر نقطه از 
نیم‌فضا و صفحه تعیین خواهند شد. از این میان 
می‌توان توزیع فشار تماسی صفحه و محیط را نیز 
تعیین نمود. هم‌چنین در حالتی که سختی صفحه 
به‌سمت بی‌نهایت میل کند (صفحه‌ی صلب) با تعیین 


نیروی کل تماسی از رابطه‌ی زیر: 


27 


۷ 2 
(۰) ب01 2272 (ه)ر۳۵ |27 ۳۲ 
ك 0 


سختی بی‌بعد فنر متمرکز معادل نیم‌فضا (شکل )٩‏ 
از رابطه‌ی زیر تعیین می‌شود: 
۲ 


3 )«( 
۷ )0(2 ۵۵ 


عبارت 2۸66 در رابطه‌ی (1۱) برای بی‌بعد شدن 
مقدار سختی اضافه له اه هم‌چنین سختی فنر 
کتتر وی معادل نیم‌فضابا استفاده از رابطه‌ی 


۵۶ 1 < ,1 به‌دست می‌آید. 


نتایج عددی 


در این بخش. نتایج به‌دست آمده از بخش‌های قبلی 


نشریه‌ی مهندسی عمرآن فردوسی 


به‌صورت عددی برآورد می‌شوند. برآورد عددی 
انتگرال‌های ارائه‌شده در بخش‌های قبلی نیاز به دقت 
خاصی دارد. چرا که حد بالای انتگرال تابم گرین 
مربوط به نیم‌فضا محدود نیست و باید با دقت مناسب 
تعیین شود. برای برآورد عددی, ۵ نوع ماده مطابق 
جدول (۱) در نظر گرفته می‌شود. در این جدول دو 
ماده‌ی اول ایزوتروپ با ضرایب پواسون ۰/۲۵ و ۰/۳۳ 
است و مابقی با رفتار ایزوتروپ جانبی می‌باشد. ماده‌ی 
ایزوتروپ با ضریب پواسون ۰/۲۵ (1 ۷12) در این 
جدول» همان ماده‌ای است که یک و همکاران [26] در 
نظر گرفته‌اند و ماده‌ی ایزوتروپ با ضصریب 
پواسون ۰/۳۳ (2 66 ماده‌ای است که لوکو و میتا 

[23] در نظر گرفته‌اند. هم‌چنین ماده‌ی ایزوتروپ 
جانبی شماره‌ی ۳ این مطالعه (3 0۷68۶ ماده‌ای است که 
اسکندری قادی و همکاران [14] در نظر گرفته‌اند. 
مطابق تحقیقات آنها سختی بی‌بعد فنر قائم معادل 
نیم فضا متشکل از ماده‌ی شماره‌ی ۱ تحت شالوده‌ی 
صلب دایره‌ای برابر ۵/۳۳۳۳ و برای ماده‌ی شماره‌ی ۲ 
برابر 7/۰۰۰۰ و برای ماده‌ی شماره‌ی ۳ برابر ۹/6۹4۸ 
می‌باشد. اندازه‌ی مناسب المان‌ها با استفاده از رابطه‌ی 
(۷) تعیین می‌شود. در این مطالعه برای صفحه‌ی 
ضلب دایره‌ای. ه برایر ۲ و ترا ضفحهی انخطاف بد یر 
برابر ۱ در نظر گرفته می‌شود. لازم به توضیح است که 
تکینگی تنش تماسی در لبه‌ی صفحه در حالت صلب 
بسیار قابل‌توجه است و نیاز به ریزتر شدن المان نسبت 


به حالت انعطاف‌پذیر دارد. بسرای تعیسین اندازه‌ی 


شنال بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 
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مناسب. با تغییر تعداد المان‌ها سختی بی‌بعد فنر قائم 
(ییب1) برای ماده‌ی شماره‌ی ۰۱ ۲ و ۳ به‌دست می‌آید 
و از آن تعداد المان‌های مناسب چنان تعیین می‌شود که 
جواب‌های به‌دست آمده منطبق بر جواب‌های موجود 
در مقالات مرتبط باشند. سختی بی‌بعد فنر قائم به‌دست 
آمده از روش ارائه‌شده در این مقاله تحت اثر شالوده‌ی 
صلب دایره‌ای با تعداد المان‌های معادل ۸۵ ۱۰ .۰ ۳۰ 
المان برای مواد شماره‌ی ۰۱ ۲ و ۳به‌دست آمده و در 
تلا ار رهق ین که 


مشاهده می‌شود در ماده‌ی شماره‌ی ۱ برای ۵المان میزان 
خطای نسبیی (مسهسکگ ‏ ( رصم 16 ۳ برابر 
۸ درصد است. در ادامه. کلیه‌ی محاسبات با در 
گفته شد. روش ارائه‌شده برای تعیین اندازه‌ی المان 
شرایط لازم برای رسیدن به دقت مناسب و در نظر 


گرفتن تکینگی فشاو تمامنی را در لبه‌های صفحه دارا 


فاد 


جدول ۱ مشخصات مکانیکی مصالح انتخابی 
5 4 3 2 1 
ایزوتروپ جانبی | ایزوتروپ جانبی | ایزوتروپ جانبی | ایزوتروپ | ایزوتروپ آمنهزه/۷ 

130000 10000 0000 5333 50000 ( تسه / ۵۲ ظ 
50000 0000 130000 533 50000 (کهه/ ۵۱ ظ 
20000 40000 20000 20000 20000 (تسدر ی 6 
20000 20000 20000 20000 20000 رکه ۵۱ 6 
025 025 025 1/3 025 0 

۱8 025 1/3 025 025 025 


جدول ۲ مقایسه‌ی سختی فائم شالوده‌ی صلب دایره‌ای ناشی از نتایج تحلیلی یک و همکاران [26] و نتایج عددی این مطالعه برحسب تعداد 


المان (ما 
اختلاف / | (60201) ,1 
۰100«( 3 > 
۱0۳۱ 3 > 
0262 0 3 > 
00112 3 > 
0037 0 3 > 
0.0015 3 > 


جدول ۳ مقایسه‌ی سختی فائم شالوده‌ی صلب دایره‌ای ناشی از نتایج تحلیلی لوکو و میتا [23] و نتایج عددی این مطالعه برحسب تعداد 


المان (ما 
اختلاف / | (60201) ,1 
5117 0 60000 
0833 0 60000 
0.0267 60000 
00117 60000 
0033 0 60000 
0.0017 60000 


سا بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


ده‌ی شماره‌ی ۱) 
(صباه) رک تعداد المان 
۰-3061 5 
۰-389 10 
۰339 15 
۰3327 20 
1 25 
۰3332 30 


ده‌ی شماره‌ی ۲) 
(«ی | تعدد امن 
۰-393 5 
ارزو 10 
۰.3914 15 
۰9993 20 
5۰9998 23 
۰39999 30 


نشریه‌ی مهندسی عمرآن فردوسی 


۱۹ بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پذیر .. 


جدول ۶ مقایسه‌ی سختی قائم شالوده‌ی صلب دایره‌ای ناشی از نتایج تحلیلی اسکندری قادی و همکاران [14] و نتایج عددی این مطالعه 


برحسب تعداد المان (ماده‌ی شماره‌ی ۳) 


اختلاف / | (60801) ,1 (حصناه) رک تعداد المان 
0230917 9.4936 164« 9 5 
02932 94936 969 9 10 
02093 94936 91923 15 
020105 94936 92938 20 
0031 94936 91945 25 
00010 94968 94947 30 
3 
اوتتا 58 ,60.21 ]۳1 هب (۲) , 
طامعصیامم 13110 صمتایتامو ۱0۲16۵1 ]سب 3 
۸/۱20 7 25 و 
۳ 2 
15 
1 
05 
1 1 0 
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شکل ۸ مقایسه‌ی تنش تماسی حاصل از حل عددی مطالعه‌ی حاضر با حل تحلیلی یک و همکاران [26] برای صفحه‌ی صلب دایره‌ای 
تحت بار متمرکز مرکزی برای حالتی که تعداد المان‌ها ۳۰ است و کوچک‌ترین المان در لبه دارای پهنای 0.0022 می‌باشد (ماده‌ی شماره‌ی ۱) 


۲ 30 
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شکل ٩‏ تنش تماسی صفحه‌ی انعطاف‌پذیر دایره‌ای در حالت چسبندگی کامل صفحه با نیم‌فضا تحت بار متمرکز مرکزی (ماده‌ی شماره‌ی ۱) 


نشریه‌ی مهندسی عمران فردوسی سال بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 
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شکل ۱۰ تنش تماسی صفحه‌ی انعطاف‌پذیر دایرهای در حالت چسبندگی کامل صفحه با نیم‌فضا تحت بار متمرکز مرکزی 
(ماده‌ی شماره‌ی ۳) 
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شکل ۱۱ تغییرمکان قائم نیم‌فضا و صفحه‌ی انعطاف‌پذیر دایره‌ای در حالت چسبندگی کامل صفحه با نیم‌فضا تحت بارمتمرکز مرکزی 
(ماده‌ی شماره‌ی ۳) 


با استفاده از حل تحلیلی شالوده‌ی صلب دایره‌ای 
توسط پک و همکاران [26]؛ تنش تماسی صفحه‌ی 
صسلب دایسره‌ای بسه‌صسورت رابطسهی 
[ 05-۳ - بت می‌باشد. همان‌طور که در 
شکل (۸) مشاهده می‌شود. روش عددی این مطالعه 
تمامی ویژگی‌های حل تحلیلی از جمله تکینگی فشار 
تماسی در لبه صفحه را به حوبی ارضا می‌کند. حتی با 
در نظر گرفتن ۵ المان رینگی برای ناحیه‌ی تماس, تنها 
حدود ۰/۵۱۰۰ درصد خطا نسبت به حل تحلیلی آن 
وجود دارد (جدول ۲). این شکل برای حالتی که 
سختی صفحه بی‌نهایت است و به‌صورت صلب 


می‌باشد. تحت اثر تغییرمکان یکنواخت ۸ ترسیم شده 


سا بیست و هفتم» شماره‌ی یک ۱۳۹۶ 


و با نمودار دقیق آن مقایسه شده‌است. همپوشانی کامل 
این دو نمودار نشان از دقت بالای محاسبات با این 
روش عددی دارد. 

در گراف‌های ارائه شده برای اهداف مقایسه‌ای, 
سختی صفحه به‌صورت نسبتی از مدول برشی آن 
نوشته شده‌است. بدین منظور مدول برشی صفحه برابر 
۴ 2 ,6 نیوتن بر مترمربع و ضریب پواسون آن 
5 2 و۷ در نظرگرفته می‌شود. هم‌چنین برای 
رعایت اختصار از این پس تغییرمکان مرکز صفحه‌ی 
دایره‌ای ((۳:)0) با ۸ نشان داده می‌شود. 

شکل‌های )٩(‏ و (۱۰) تش تماسی صفحه‌ی 
انعطاف‌پذیر دایره‌ای و نیم‌فضا را با در نظر گرفتن 


نشریه‌ی مهندسی عمرآن فردوسی 


چسبندگی کامل صفحه و نیم‌فضا (<0) تحت 
نیروی متمرکز با اندازه‌ی واحد در مرکز صفحه برای 
مواد شماره‌ی ۱ و ۳ نشان می‌دهد. همان‌طور که در این 
گراف‌ها مشاهده می‌شسود تنش تماسی برای مقادیر 
مختلفی از سختی صفحه رسم شده‌است و طبق انتظار 
پا کاهش سختی صفحه. تنش در مرکز صفحه زیاد 
می‌شود و با افزایش سختی آن مقدار تنش در لبه‌ی 
صفحه زیاد می‌گردد. شکل‌های نشان داده‌شده نشان از 
حاکم بودن رابطه‌ی (۵۲) در حالت تماس کامل صفحه 
و نیم‌فضا را دارد. شکل (۱۱) تغییرمکان قائم نیم‌فضا و 
صفحه‌ی انعطاف پذیر دایره‌ای را برای ماده‌ی شماره‌ی 
رای هاگن کاس مس وتف 
سختی‌های مختلف صفحه تحت بار متمرکز در مرکز 
صفحه نشان می‌دهد. با توجه به این شکل مشاهده 
می‌شود. با کاهش سختی صفحه تغییرمکان انتهای 
صفحه بیشتر می‌شود و با افزایش سختی آن, مقدار 
تغییرمکان صفحه به تغییرمکان مرکز آن نزدیک‌تر 
می‌گردد. هم‌چنین به‌دلیل برقراری چسبندگی کامل 
تک وی قضا ظی رای (۵1) اشتان یرتک 
قائم صفحه و نیم‌فضا در کل شعاع صفحه با هم برابر 
شا تنل 

همان‌طور که پیشتر گفته شد. مقدار شعاع 
که نشج طوآن نک میرن اف کوش 


1۱23 
۳0 1۴ 0 

1 / )6, 0*( 2 1 

۸۱-0 


بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پذیر .. 


روند تحلیل تعیین شود. نحوه‌ی تعیین آن طبق شکل 
شماره‌ی (۱۲) به این صورت است که هم‌زمان با 
کاهش مرحله‌ای شعاع صفحه از 2 تا حدود 0. مقدار 
تتش تماسی در هر مرحله ترسیم می‌شود. این روند تا 
یافتن اولین شعاع صفحه که به‌ازای آن» تنش در تمام 
نقاط صفحه از شعاع صفر تا آن شعاع (0) همگی 
مثبت باشند. ادامه می‌یابد. در جدول (۵) مقادیر 0 
به‌دست آمده از این مطالعه با مقادیر به‌دست آمده از 
نتایج یک و همکاران [26] برای ماده‌ی ایزوتروپ 
شماره‌ی ۱ به‌ازای سختی‌های مختلف صفحه‌ی 
انعطاف پذیر دایره‌ای» مقایسه شده است. اختلاف ناچیز 
این دو مطالعه نشان‌دهنده‌ی دقت بالای محاسبات 
می‌باشد. هم‌چنین در جدول (1) به‌ازای سختی 
1 2 (*0,7)/ 1 برای صفحه‌ی دایره‌ای. مقادیر .0 
ب‌دست آمده برای مواد ایزوتروپ و ایزوتروپ جانبی 
با هم مقایسه شده‌است. با توجه به این جدول. مقدار 
0 در شرایط یکسان برای ماده‌ی شماره‌ی ۳ کمتر از 
ساير مواد می‌باشد و بدین معناست که ناحیه‌ی بلند 
شنلده‌ی ضفحه در این ماده پیشتر از سایر مواد اس که 
نشان‌دهنده‌ی سخت‌تر بودن این ماده است. از طرف 
دیگر ماده‌ی شماره‌ی ۲ دارای بیشترین نرمی بین مواد 
قوش 
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ن ۲ روند یافتن شعا حسیند 5 صفحه و نیم‌فضا (0 برای صفحه‌ی انعطاف یذ دار ه‌ای در حالت عد حسیند 5 صفحه با نیم‌فضا 
تن سعاع ۰ سم بر بر ابر 231 تیم 


تحت بار متمرکز مرکزی (ماده‌ی شماره‌ی ۳) 


نشریه‌ی مهندسی عمرآن فردوسی 


شنال بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


۱۹ 


جدول ۵ مقایسه‌ی شعاع تماس (0) به‌دست آمده از این مطالعه با نتایج پک و همکاران [26] به‌ازای مقادیر مختلف سختی صفحه برای 


صفحه‌ی دایره‌ای انعطاف‌پذیر (ماده‌ی شماره‌ی ۱) 


اختلاف 1 | [۵1]2008 6 2۵ مطالعه حاضر ((6,۵)/ظ 
20343 0870 09572 005 
25 01590 033 003 
22508 06220 06080 0015 
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جدول 1 شعاع تماس (0۸) به‌دست آمده برای مواد مختلف نیم‌فضا به‌ازای سختی 0.01 2 (۵7 ,1(/)0 صفحه‌ی دایره‌ای انعطافپذیر 
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شکل ۱۳ مقایسه‌ی تنش تماسی صفحه‌ی انعطاف‌پذیر دایره‌ای با نیم‌فضا در دو حالت چسبنده و عدم چسبندگی صفحه با نیم‌فضا تحت بار 


متمرکز مرکزی (ماده‌ی شماره‌ی ۳) 
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شکل ۱۶ مقایسه‌ی تنش تماسی صفحه‌ی انعطاف‌پذیر دایره‌ای با نیم‌فضا بین مواد مختلف نیم‌فضا در حالت عدم چسبندگی صفحه با نیم‌فضا 


تحت بار متمرکز مرکزی 


سال بیست و هفتم شماره‌ی یک ۱۳۹۶ نشریه‌ی مهندسی عمران فردوسی 
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بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف پذدیر 5 


شکل ۱۵ مقایسه‌ی تغییرمکان قائم صفحه‌ی انعطاف‌پذیر دایره‌ای و نیم‌فضا در دو حالت چسبنده و عدم چسبندگی صفحه با نیم‌فضا تحت 


بار متمرکز مرکزی (ماده‌ی شماره‌ی ۳) 


0 قیازیسی رون تهالا-ق 
چسبندگی کامل صفحه‌ی دایره‌ای با نیم‌فضا و هم‌چنین 
عدم چسبندگی آنها برای ماده‌ی شماره‌ی ۳ و سختی 
ثابت 0.01 (*9/)0,۵ ترسیم شده‌است. هم‌چنین 
در شکل (۱۶) تنش تماسی در حالت عدم چسبندگی 
صفحه‌ی دایره‌ای با نیم‌فضا با همان سختی صفحه بین 
مواد مختلف نیم‌فضا به‌صورت مقایسه‌ای ترسیم شده 
است. اثر میزان ناهمسانی نیم‌فضا بر پاسخ نیروی مژثر 
بر صفحه‌ی انعطاف‌پذیر مشاهده می‌شود. شکل (۱۵) 
تغییرمکان قائم صفحه‌ی دایره‌ای و نیم‌فضا را در دو 
حالت چسبندگی صفحه با نیم‌فضا و عدم چسبندگی 
آنها با مقدار سختی ابت صفحه نشان می‌دهد. 
همان‌طور که مشاهده می‌شود در شرایط عدم 
چسبندگی, از شعاع صفر تا 0 صفحه و نیم‌فضا با هم 
در تماس‌اند و از شعاع 0 تا ه صفحه از نیم‌فضا جدا 
خواهد شد و تغییرمکان قائم نیم‌فضا از صفحه بیشتر 
می‌شود که نشان‌دهنده‌ی حاکم بودن رابطه‌ی (۵۳) 
می‌باشد. 

شکل (۱۱) مقادیر تتش قائم نیم‌فضا را در سطح 
نیم فضا به‌ازای مقدار سختی ثابت صفحه دایره‌ای نشان 
می‌دهد. همان‌طور که مشاهده می‌شود. با دور شدن از 
مرکز صفحه. طبق اصل منظم بودن جواب. مقدار تتش 
قائم نیم‌فضا به سمت صفر میل می‌کند. شکل (۱۷) 


نشریه‌ی مهندسی عمرآن فردوسی 


مقادیر تغییرمکان قائم نیم‌فضا در عمق را به‌ازای مقدار 
سختی ثابت صفحه‌ی دایره‌ای نشان می‌دهد که در این 
مورد نیز مطابق فوق با دور شدن از محل اثر با مقدار 
تغییرمکان قائم صفحه به سمت صفر میل می‌کند. شکل 
(۱۸) تأثیر میزان سختی صفحه بر تغییرمکان قائم 
نیم فضا در سطح را نشان می‌دهد. این شکل در حالت 
چسبندگی کامل صفحه‌ی دایره‌ای با نیم‌فضا به‌ازای 
سختی‌های مختلف صفحه ترسیم شده‌است. همان‌طضور 
ها متسر توف اف انم تس ,حون 
یکنواختی تغییرمکان قائم نیم‌فضا در سطح زیرصفحه 
| 
تغییرمکان غیریکنواخت و روند کاهنده خواهد داشت. 
در انتها لازم به توضیح است که با توجه به 
کات هاش ناوات تاه ان تا نت نو تروش قو 
پاسخ مسأله بدین صورت است که مقدار شعاع تماس 
() بانسبت ,0/ وو۸ رابطه‌ی معکوس دارد 
به‌ طوری‌که با افزایش این نسبت. شعاع چسبندگی 
صفحه کاهش می‌یابد. هم‌چنین روابط ارائه‌شده در اين 
مقاله. برای نیم‌فضای ایزوتروپ جانبی می‌باشد و با 
انتخاب مصالح ایزوتروپ نمی‌توان از این روابط 
استفاده نمود چراکه برای نیم فضای ایزوتروپ نیاز به 
بازنویسی روابط و اصلاح مقادیر ,5 و ره و ساده‌سازی 
روابط می‌باشد که در مرجع [26] ارائه شده است. 


شنال بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 
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شکل ۱۰ مقایسه‌ی تنش قائم نیم‌فضا در سطح بین مواد مختلف نیم‌فضا در حالت عدم ج چسبندگی صفحه‌ی انعطاف‌پذیر دایره‌ای با نیم‌فضا 
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شکل ۱۷ مقایسه‌ی تغییرمکان قائم نیم‌فضا در عمق بین مواد مختلف نیم‌فضا در حالت عدم چسبندگی صفحه‌ی انعطاف‌پذیر دایره‌ای با 
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شکل ۱۸ مقایسه‌ی تغییرمکان قائم نیم‌فضا در سطح به‌ازای سختی‌های مختلف صفحه در حالت چسبندگی صفحه‌ی انعطاف‌پذیر دایره‌ای با 


نیم‌فضا تحت بار متمرکز مرکزی (ماده‌ی شماره‌ی ۳( 


سال بیست و هفتم شماره‌ی یک ۱۳۹۶ نشریه‌ی مهندسی عمران فردوسی 


۱6۲ 


نتیجه گیری 

در اين مقاله» به تحلیل توأّم نیم‌فضای ایزوتروپ جانبی 
و صفحه‌ی دایره‌ای انعطاف‌پذیر تحت اثر نیروی قائم 
متقارن نسبت به محور گذرنده از مرکز صفحه پرداخته 
شد. برای انجام تحلیل ابتدا معادله‌ی دیفرانسیل حاکم 
بر صفحه حل شد و از آن تابع گرین تغییرمکان قائم 
صفحه تحت بار حلقوی وارد بر آن به‌دست آمد. سپس 
نیم‌فضای ایزوتروپ جانبی تحت اثر بار حلقعوی فائم 
وارد بر آن بررسی شد و با حل معادله تعادل حاکم بر 
محیط. تابع گرین تغییرمکان قائم نیم‌فضا نیز به‌دست 
آمد. در ادامه با در نظر گرفتن شرایط پیوستگی صفحه 
با نیم‌فضاء تابع تنش تماسی صفحه از حل معادلات 
انتگرالی حاکم بر مسأله با استفاده از روش المان‌های 
رینگی تعیین شد. با داشتن این تابع و انجام مراحل 
سعی و خطاء شعاع چسبندگی صفحه با نیم‌فضا تعیین 
شد و از نتایج آن نمودارهای تنش و تغییرمکان برای 
صفحه و نیم‌فضا در سطح تماس و در نواحی 
دوردست به‌صورت گراف‌هایی ارائه شد. نمودارهای 
ارائه‌شده نشان داد که با کاهش سختی صفحه‌ی 
انعطاف پذیر مقدار تنش تماسی (0,)۲ در زیر محل 
اثر بارگذاری روند افزایشی دارد به‌طوری که هرچه 
سختی صفحه به صفر نزدیک شود مقدار تنش تماسی 
در محل بارگذاری به سمت بی‌نهایت میل می‌کند. از 
طرفی دیگر با افزايش سختی صفحه مقدار تنش 
تماسی در لبه صفحه بیشتر می‌شود به‌طوری‌که در 
صفحه با سختی بی‌نهایت (صفحه صلب»» تنش در لبه 
آن به‌سمت بی‌نهایت میل می‌کند. هم‌چنین مشاهده شد 
که اثر میزان ناایزتروپی در پاسخ مسأله به نسبت 
6 / و۸ وابسته است به‌طوری که با افزایش نسبت 
فوق. جداشدگی صفحه و نیم‌فضا در شعاع کوچک‌تری 
اتفاق می‌افتد. 


مرا 
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نمودارهای مربوط به تغییرمکان صفحه‌ی 
انعطاف پذیر و نیم‌فضا نشان می‌دهند که مقادیر آنها 
تابعی از جنس مصالح صفحه و نیم‌فضا می‌باشد و اثر 
متقابل سختی این دو محیط. بر آنها تأثیرگذار خواهد 
بود. به‌طوری که با افزایش سختی صفحه مقدار 
تغییرمکان صفحه و نیم‌فضا یکنواخت‌تر می‌شود. با 
مقایسه‌ی اثر مواد مختلف نیم‌فضابر تش‌هاو 
تشیترشکان ها نخان داده ند کدی مانهی ۳۲ سفن 
پیشتری نسبت به ساير مواد دارد و بالعکس ماده‌ی 
شماره‌ی ۲ سختی کمتری نسبت به ساير مواد دارد. 
هم‌چنین با توجه به جداول و گراف‌های ارائه‌شده 
مربوط به محاسبه‌ی تنش تماسیء ابت شد که با 
کاهش سختی صفحه‌ی انعطاف‌پذیر شعاع تماس 
صفحه و نیم‌فضا (0) در نقاط نزدیک‌تری نسبت به 
محل اثر بارگذاری اتفاق می‌افتد و با افزايش سختی 
این شعاع به لبه‌ی صفحه نزدیک‌تر می‌شود به‌طوری که 
در حالت سختی بی‌نهایت. صفحه هیچ‌گاه دجار 
بلندشدگی نخواهد شد. هم‌چنین پدیده‌ی جداشدگی 
در گراف‌های مربوط نشان داده شد که طبق انتظار از 
شعاع بعد از 0 تغییرمکان نیم‌فضا از صفحه بیشتر شد 
و با پلند شدن صفحه عدم پیوستگی این دو محیط در 
این ناحیه مشاهده شد. در انتها با توجه به شکل‌های 
تغییرات تنش و تغییرمکان نیم‌فضا نسبت به سطح و 
عمق در نقاط دوردست. مشاهده شد که با دور شدن از 
محل اثر بارگذاری مقادیر آنها به سمت صفر میل 
می‌کند که دلالت بر اصل منظم بودن جواب دارد. یکی 
دیگر از نتایج قابل ذکر در این قسمت. محاسبه‌ی 
شش تفایکن شعاد عیفر خای است 
که سختی صفحه بی‌نهایت است (صفحه‌ی صلب) که 
با مقادیر تحلیلی موجود به‌ازای مواد مختلف نیم‌فضاء 


مطابقت دارد. 


رحیمیان محمد و اسکندری قادی» مرتضی ‏ "تئوری ارتجاعی " انتشارات دانشگاه تهران. (۱۳۷۹). 


۲ . اردشیر بهرستاقی. عزیزاله و اسکندری قادی» مرتضی . "تحلیل نیم‌فضای ایزوتروپ جانبی تحت اثر صفحه‌ی صلب دایره‌ای با استفاده 


نشریه‌ی مهندسی عمرآن فردوسی 


شنال بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


قاسم گرجی بندپی- مرنضی اسکندری قادی- عزیزالله اردشیر بهرستاقی 


از توابع گرین بار حلقوی" نشریه‌ی علوم کاربردی و محاسباتی در مکانیک, ۲۲ (۱)؛ صفحات ص.ص. ۵۸-4۳ (۲۰۱۱). 
:1 ,(114)2 ,.1۷66۳ .۳ .ل و"20108صا0] متاقواه ووماصمعمعه مه عصااوع۲ ومفقام تقانم‌ممامی؟ ,بر روعامن 
(1988) ,2083-2092 00۰ ,([14:12)2083 1/)۵۸0۳(0733-9399-1 10۰106 
:1 ,(114)10 .۷۲66 .9 بل رصملاه0صام؟ تمللصد/۳ وعمامماعمه؟ صم متقام تقلتمتت .2 وان 
.(1988) ,1723-1739 .00 ,(0733-9399-114:10)1723 ]10۰1061/۸8 
وععلصمتعمه 2 جم صصههها وه ۵۶ مممتاه:۷۱۵ ۳۳۵۲۵۵0 وی رصلهعوت تاه 4ص2 ربق مقعلله۷21( و2 محعامن 
,(1989) ,235-246 .00 و2 .15516 و128 ,۷۵۱ ۷۱۵۳۵۵۵0 ۵80 500 ]۵ آمد ناه ۲ +صمتاه‌صنام] 
۵ 509۵۱0 ]۵0 09۲۳۵۱ .۲۵۱0021109 ووماطمتعجه) ج م۵ ففصید تقانمتن ۵۶ فممتاه:۷10 مصقام-ص1 رد2 رحعامن 
1990(۰) ,461-471 .00 و3 .5506 ,143 ۷۵ ر1۷۵۲۵0۳۲ 
0۲ متاوهاه تمامجصعهم-7۵ ووماممتعصم) 2 جم حصهه‌ها م۲۱ تهلناه‌ت۳ ۳.۰ وتتصعنا 20 ور رحواعن 
.(2005) ,431-439 .0ج و6 .55116] و85 ,۷۵۱ رکعز۷۵6۳۵۳ 00 ۱۷۵۵۱۵۱۵۵ ۱۵0 آمم ]۵ امسو7 ]7۸0 
۶0۱0۱201010 وومطم1قصهعا ماممصونه۳۷۵0-۵ 2 جم صنهوها ۱۵2060 والهمتمصتر" ۲۰ رتتصونا 224 ور بقعاعن 
.(2007) ,555-574 .00 و5 0۰ و27 ۷۵۱ روعتصهطمع]۱۷ 0ج مصتنممصموص امتتامنانه 
2 عم ماهام تقایمتته وعطظ متاعقاه مه ۵ مصماهط ۵0ع10۲ ماتامصصططترونظ نک رتقانات مضه وم روعامن 
۰ ,3-5 .1851165 و301 ۷۵ ۷۱۵۲۵۵۵ ۵0 یمه ۵ هیال وط0اجصا0؟ تمامحوتوم ۵ وو1۵ظ810ها 


495-509, )2007(. 


7۲ ملظ .عاهتتمممط مامرمتاموا وامقه7عصهط که وااملافهاط ربا رعصقطا2 فصه ,۷۷ رحعطن 20 یط رعطان . 


2006(۰) و۵۲[205ظ)۱۵ عظ 1 10۲0۲60 ۸۵ 3300 و17 :5۲0 ۲.۵۰ ر6ع80110 


و"کلقاوه 16حمتاماممه اقطمعمعمط ما مملاتاطتتاونل . ققمه . اقصمتعممصصاه.. معط .رنه بل لاملا . 


(1948) ,522-533 .00 ,44 ,۷۵۱ ,۹001607 اهمنطممومانطط معلتتطاصصون مطا ۵۶ ومصتلععع۳۵ظ 1همتاممومط)۱۷]2 


,"۵12 1901070010 «اموه7قصج 10۲ فممتاهنا60 ۷۷۵۷۵ مطا ۵۶ فممتتاامو مامامصمع ۲۸ ۷۲۰ رتهدان-021صق او 


2005(۰) و1-19 00۰ و81 ,اکتا ]0 .ل 


000 ۲2۱6۷۵۲۹۵17 ]۳ ریم مبتطمهاگمتطهظ دوع مضه ویو ,۷۲ نک بلط و بتقطن21-0صم او . 


5 901 ,"10205 عمهناو تملمصا اناد ماتصصا جح وم تمع1 متص ه ۶ مهو ملصصح‌صرول‌ماکهآه 


.(2008) ,986-1003 .00 و8 عصاتهمصلوص۲ معلمنمطنه 290 


0 0۴ ممتاهط۱1 ۷۵۲62 ۲۳۵۲۵۵0 بیط رتطعه)ومتطه نومه مه و۷ رتطعاله۲ و۷ رتلقطان- لصو .. 


,۷ لک ۵ ۱۱6۵۲۵6۵ 2۳91۱6۵۲۱9 ]0 21-502000۲۵ مامماوو1 راوو7۵فصجتا ه وم مفلل تهلتام‌تام 


00. 913-922, )2010(, 


م مامت منود ۵۶ ممتاهه0 اممتا۷۵ ۳۳۵۲۵۵0 ریم رلطعهاومتطه عنم را تقفطان-هه ماو .. 


ماه 40ج ومتصصمطو ام موی کاقط مدممتا0و1 ام‌ونهتعصه چ ۵۶ طاممل بحتهتازمانج مج صا 4عتتناه 


2010(۰ ,547-560 .00 ,ز/10.1016 :مک ,30 عصنمهع‌صنعوه 


40 ده ۵۴ ۷۱۵۳۵2610 مصلاهم؟۴ به راومه طوم0ته مضه رم رتتا0م۱۷۱۲22 ری تقطن021-0صم او .. 


0 مامح اهعت) قمع صرح لهمزمصصی رممجوی لا مزممتا0و1 آم‌وه7عصه چ صه صا ۵0زناه وونل تقلتام‌تام 


(2010) ,547-560 .0 ,976 .10۰1002/88 :401 روعتصهطهمعطط0وع 


همهم رلقصح ۱/۵۲۱۵ ۳۸ رم رتطعماوم‌تطهظ- ددع۸0 200 ریگ ۷۰ بظ ئلهظ و۷ لصو . 


هط متاعهاه 1501۲0۵16 موه عصها ج ره عمجم طمصتی تهاناصها نهر متاجممومم لهع۷۵۲ ۶۵۲ 2۵00 


(2011) ,10.1115/1.4005546 :12۵۲ 1۸۱۷۲ :۲1 _ژ رکع۱۱6۵۳۵۴ ۱60 هام0 ]۵ 0۱۲۵1 ,502606۳ 


۱5۳ 


سال بیست و هفتم شماره‌ی یک ۱۳۹۶ نشریه‌ی مهندسی عمران فردوسی 


۳32 بررسی تحلیلی اند رکنش صفحه‌ی دایره‌ای انعطاف‌پذیر .. 


یلع ه ۶۵۲ وتمعطا راتمتافهاه ۵۶ مصمامامنص متتامصصتنرفلنه فص مر ری راما 24 ربق نک رقکلصوطانا .. 
1954(۰) ,89-101 .00 رکخکاموه 00 م۱۱66 آمم۲۵8۵ 0۲ ما0 150000 ۲7۵158 فصعتا ۱۳ 

لفط متاوجاه مج 0ص ماقم تماتامته ۵ ممع۲اعه امماصمن ۵0۵ص ویکگ. یک وتعت 4صه ربآ0 01207 
1974(۰) ,115-130 00۰ و4 ,۷۵۱ ,56اب ]0 0۱۲۵1 ,"5026۵ 

,۷12/2 ۷۵۱۰ ۳۳5۱ تم طمنامال‌صه۲۱ و(:۵0) معمیاا۳ .9 نصا و واتمتافهاه ۵۶ رمقط) تحعصا فظ ۳ ر.ط مصتاتتت 
.(1972) ,1-295 00۰ ,۷0116 ۲6۲ ۲۱۵۱00106۲۵ صناتوظ رته‌عصنمو .۲۳۵80611 ,0 ۵0 رل 90105 ۵۶ ومنصهطهع]۱۷ 
0006 اوه *عصهتا و ۵۶ واتهتاعهاه ۵۶ معط فص ۶ه عمرمامامعج آقطمتعصمصصتل معتط مطه م۵ وت ۲ وتا 
1953(۰) و145-151 .00 و2 ,۷۵۱ رهعلصلگ لو و۵0 

۳۲۵615160 صقو ریم عصتطفتامانام وجر-صومام۲۱ رد وعزلهها متاقفاه مام۲0امونصع ۵۶ ۳1۳۵۵7 ,.ت .9 رللاعاتصطام ][ 
19861(۰) ,]0211 

متاقعاه ما ممهموفط صرمنصنا چ من صمتاه‌صام؟ تقاتمتام وچ ۵۶ مقممم‌ومگ. ربق مهال۱ 220 ور .1 معا[ 
.(1987) ,105-118 .00 و15 ۷۵۱۰ و۱265 آمتامتناه 4صه موصعم ملمتمطاته 7۷2۵/۵۵ 

08 مر 0۱01277 مصقام متصصا مه ۳1 50110 متافهاه 0۲0۵10مام2 صیح صا ووعتاو م۰1" ویلط .1 مالقطم۷1( 
.(1900) ,247-258 00۰ ,32 ,۷۵۱ ,5001607 1621)محصمطممه مل‌صم؟ مطا ۵۶ 

7 متافقاه صه عومجم فصصمامامم اقصمتعصمصتل عععطا صا صمتامصیظ ققواه فا .۲۷۷ رلکامه0۳ 
.(1954) 21-25 00۰ و2 ۷۵۱۰ و56 ۳۵۱۵۲۰ .20ع۸ بلاط 1500۲00 «اووزه7عصهتا رها ۲1200ماممتهده 

لمح منقام ماطلععا؟ ه ۵۶ امفاوم ووملصمنفصهآ رت .ل رعل6ملطفه 20 ری .ظ رعصمصصصصا ریق ۷۰ بظ ملظ 
۱6۳۵۵۵۳۵ ر6ط0اجناوه لهتعهاصا رها و1۵20 مرتامصصمهررعتنج تعصتا ممهوری هط طاممصصی ج صز کنانتصصج 
(2008) ,1004-1011 .00 ,50 5668665 ۱۱6۵۳۵۳۴6۵1 ]0 0۱0۲۵1 

006 01و6۲ ۲عصها ماتصاکص‌تصمی نم فممتبامو عفممتامص قصممز و۲۷۷۰ ۲۰ رتامطن) 280 وین ۷۰ معط 
.(1979) ,545-551 .00 و(17)5 ,56 .5۳9 . 16 ,"۳0۵16۲1215 

1951(۰) ,۷ .۲ وال۷۵ ۲و الا ۷۱۵0۵۲ و عصمکعصهتا تفتتتام۳۴ بلط :] رطمل8860 

و106 ما۲۱ ۳۳۵۵۱۵6 ۱۵0008 لهمزرمصصیل صرح [2م1وعهان رقعاهان 0۶ فلوولقصه 2240 ۰۲۳۵۵۲۲" ریبک 821127 
.(1974) ,6۳50617 ۱۵ رکنان ۳8۵127۷۵00 

نعلر ۷ ۱۱۵۲ .60 2۳ ,"علامطه 4صح فعفقام ۵۶ معط رگ رهلک او مصز۷۷0 فصع ...6 ومعلص‌طومص] 
.(1959) ولد ۱۷۲۵0۲۵۷۷ 

۷ له فصمتاتامی امتمصهع ۵۶ ووعصمیمتصاطمط 4ج وومصمامآم‌مط۵ .۷۷ رعصه ۱۷۷ 200 و2 ۷۲ رعطه۲۷۷ 
.(1995) ,501-513 .0 و(374) 32 رکععت اه وه ک8نا6 ]۵ آمص سول ۲۱۵۲۳۹۵۲۵0۳0۵1 و امتاوهاه 1501100160 
۲۵اه "م2٩‏ کاقط متاعهاه جح صح عمتقهام ممع07و آمماومن 060صنان‌طاصنا ما پا ولا رطقططوازه ۲۲ 
.(1969) ,198-202 .00 ,36 :۸9۱۳ روی۱۷6۵۳۲۵۳ ۱60 امه 

2 0۲ص )میم متاقعاه فعماممتعصه ه مضه صصهعه ه ممم07وه مملامتجم۵و ۵۶ امقصم ولا رطهصصعاز ۱۷۷ 
1971(۰) 707-711 00۰ و(8) 13 ,56۱6۴665 ۱۷۵۵۳۱۵۴۵6۵1 ]0 یال 1۱6۲۳۵۵0۳۴۵1 1020 ۳2۵۷1 
6۲۵۵۵۵ ,"211-50206 متاوهاه صح مضه حصهوها ه ممع07عه عمجم ووملصمنعصها ۸ ول رصهصصیازه ۱۷۷ 


-(1972) ,00.73 و10 :۷۵1 ,56۱6۳166 0۲1۳۳9۱۳۵۵۲9 0۱۲۲۱۵۱ 
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